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ABSTRACT 

We show that accretion disks with magnetic fields in them ought to make jets provided 
that their electrical conductivity prevents slippage and there is an ambient pressure 
in their surroundings. 

We study equilibria of highly wound magnetic structures. General Energy theo- 
rems demonstrate that they form tall magnetic towers whose height grows with every 
turn at a velocity related to the circular velocity in the accretion disk. 

The pinch effect amplifies the magnetic pressures toward the axis of the towers 
whose stability is briefly considered. 

We give solutions for all twist profiles &(P) — fi(P)i and for any external pressure 
distribution p(z). The force-free currents are given by j = 5(P)B and we show that 
the constant pressure case gives a cx which leads to analytic solutions for the 

fields. 

Key words: Jets, quasars, star-formation, MHD. 



To take some of the most dynamic objects in the universe 
in which 'apparently superlummal' motions have been seen 
and to work on such problems using statics may seem an in- 
dication of a seriously deranged scientist. Nevertheless, the 
potential energies involved in any problem can be studied 
statically and it is those potential energies that drive the 
motion. Thus, without studies of the way the potential en- 
ergies operate, the basic understanding of why the motion s 
occur in the way they do may be lost. As Eddington (1926) 
said 'the chief aim of the physicist in discussing a theoretical 
problem is to obtain "insight" - to see which of the numer- 
ous factors are particularly concerned in any effect and how 
they work together to give it'. 



Even a perfect model of a phenomenon, that gives all 
the observables correctly, is not good science until it is 
analysed to show which aspects are essential for the phe- 
nomenon. 

An unnecessarily detailed model, which reproduces the phe- 
nomenon, can actually be a barrier to understanding. The 
lack of real understanding of what makes red-giants is a 
case in point! It may require the insight of an Eddington 
rather than the calculations of a Chandrasekhar to simplify 
the model to the bare essentials. 

Within an accretion disk any radial magnetic field will 
be sheared and stretched by the differential rotation, so the 



resultant toroidal magnetic field will grow until it is strong 
enough to arch up out of the disk with the gas flowing back 
down. 



Hoyle & Ireland (1960) were the first to describe this 
instability but the con ditions for its occurrence were more 
precisely calculated by Parker (1965). In the resulting con- 
figuration the mass of gas in the accretion disk continues 
to anchor and to twist the feet of the flux tube. We show 
here that if there is an ambient pressure in the tenuous gas 
above the accretion disk then the flux tube will grow in 
height with every twist of its feet but will not expand lat- 
erally. Thus a tall tower of magnetic field is formed whose 
collimation or aspect ratio Z/R increases linearly with the 
number of twists. 

My thesis is that an important pre-requisite to under- 
standing the dynamics of the jets above accretion disks is a 
serious study of the magneto-statics of the magnetic fields 
that they twist into their coronae. My first studies in this 
direction ended in total failure. In 1979 I had a mechanism 
that would give a tall tower of magnetic field that grew more 
collimated with every twist. With much enthusiasm I made 
a more detailed and exact calculation but to my amaze- 
ment and chagrin the field managed to expand to infinity 
and disconnect itself after just over half a turn! The hoped 
for collimation that got better with every turn ended with 
only half a turn when the degree of collimation was not a 
needle-like, one degree, but a full 120 degrees! See Figure 1. 
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Figure 1. After half a turn a field stretches to very large distances 
in the absence of a confining coronal pressure. After -^L turns = 
207 deg the field stretches to infinity and further turning does no 
work Fo r more recent work on field opening see Uzdcnsky ct al 
(2002Qj]nd Uzdcnsky (2002 ). 



Years later I wrote up this problem ( Lynden-Bell fc Boily 
1994 Pi.per I) for a conference celebrating Mestel's work. It 
pr oved to be a fine e xample of a phenomeno n long advocate d 
by |Aly (1984, 1995j ) in solar MHD, see also |Sturrock (199l[ ). 
This rekindled my interest, but it still took two years for 
me to realise an ambient coronal pressure external to the 
field would prevent it expanding to infinity. To do so would 
take too much energy. When the field cannot expand to in- 
finity, the continual turning of the accretion disk does wind 
many turns into the corona, provided the conductivities are 
great enough for flux freezing at the feet. Then the original 
argument comes into operation and tall towers will be gen- 



understanding of this problem for any specified fi(-P) and 
any chosen external pressure field p(z). We emphasise that 
we have the simplest magnetically dominated model. Above 
the disk, magnetism so dominates that the field is force- 
free everywhere within the towers, but at their surfaces the 
magnetic forces balance the ambient hot gas pressure p(z). 
Where there is gas there is no field, where there is magnetic 
field there is no gas pressure but we nevertheless assume 
perfect electrical conductivity. When we view our sequences 
of equilibria parameterised by t they are so reminiscent of 
what is seen in radio galaxies, quasars, star-forming disks 
with Herbig-Haro objects, etc., that it is hard to resist the 
temptation of talking in the terminology of dynamics rather 
than statics. Please remember that the velocities we talk 
of are only the velocities that would occur in the absence of 
inertia when accelerations can become arbitrarily large with- 
out penalty. We consider som e effects of inertia in paper IV 
(Lynden-Bell, in preparation). 

We believe that the very simple calculation given in sec- 
tion 2.3 contains the answer to the question 'Why do flat ac- 
cretion disks produce needle like magnetic jets?', but further 
refinements, given in paper IV on detailed field structure, 
give greater understanding of why jets act as giant linear 
accelerators and why their electric currents are so concen- 
trated to the tower's axis. 

In the above, I have given a personal account of the 
evolution of my thoughts and calculations but there are 
many other ideas in this field, which started with the Jet 



in M87 observed by Curtis (1918). Although Schott (191; 



wrote a book which detailed the theory of synchrot ron ra- 
diations in 1912, it was not applied to M87 until |Baadc 



erated whose heights grow with each successive turn. Basic (1956) ). Even though the first quasar 3C273 had a promi 



t heore ms on magneto-statics and much simplified models of 



such m agnetic towers in a constant external pressure were 



given in paper II (Lynden-Bell, 1996) which used a rigidly 



nent j et these remained an enigma emphasised by Wheeler 
(197l|) in 1970, who ri g htly drew at tention to a paper by|Le 



rotating inner disk and a fixed outer disk to which the field 
ret urne d. Simplified models with realistically rotating accrc- 



Blanc fc Wilson (197Cj ); IRees ( 1971ft provided one of the firs t 
expl anations which evolved into Plandford fc Rees (1971) 



tion disks were derived in a conference paper Lynden-Bell 
(20010 The present paper is a development of those calcu- 
lations to allow for a pressure that decreases with height. 
We also find a better approximation for the distribution of 
twist with height along each field line. This improves on the 
assumption of a uniform distribution made earlier. 

We derive the magneto-statics of force-free magnetic 
fields whose feet have been twisted by an accretion disk. We 



and jcheuer (1974 ). Barnes & Sturrock (1972) and Lovelace 
(1976| ) were among the first to t reat magnetically dominated 



models, while for black holes Blandford & Znajok (1977 
gave such ideas a new and most interesting twist by produc- 
ing a magnetic me chanism for extracting the spin energy 
from a black hole. Blandford & Payne (1982) have a nice 



mechanism for driving winds centrifugally. Extra-galactic 
radio jets were reviewed with a wo n derful series of radio 
pictures by Begelmann et al. (1984). Sakarai (1985, 1987) 



showed a slow asymptotic collimation of s uch winds and the 



magnet 


ic field lines labelled P being twisted by an angle 


wind equations have been much studied (Heyvaerts & Nor- 


*(P) = 


Q(P)t at their feet which are anchored in the disk. 


man, 1989; Appl & Camenzind, 1993; Pudritz & Norman, 


Here Q 


P)t is the differential rotation of the feet of the flux 


1987; |Lovelace et al., 1991). pkamoto (1999) criticised the 



tube and f is a parameter that gives the amount of that 
twist. If we increase t we change our equilibrium along a 
Poincare sequence. Thus, in the language of his catastrophe 
theory, t is our control parameter. Of course, if all inertia 
were unimportant and if a real accretion disk were drag- 



collimation claims of most work ers. Computational simula - 
tion s of jets have been m a de by [ Bell fc Lucek (1995, 1996 ) 



Ouycd et al. (1997j ); puyed fc Pudritz (1999| ). Many 



ging th 3 feet of the flux tubes at relative angular velocity 



f2(P), then this sequence of equilibria would give us a film 
of how the field would evolve in the presence of the external 
pressure field p(z). This would not give a true picture just 
because the inertial terms are not always negligible, nev- 
ertheless such sequences are very instructive in that they 
show us how the field would like to change if inertia did not 
slow its accelerations. We find it possible to get a reasonable 



and 

of these stu dies start with a unif orm magnetic field at in- 
finity as did Lovelace (1976] ) and Shibata fc Uchida (1985^ 
1986j ) who produced very convincing jets this way. However, 
a uniform field at infinity puts in a collimation at the start, 
whereas some wish to see both the field and the collimation 
to emerge as a consequence of the persistent rotation rather 
than having them inserted as a boundary condition. Jets 
are by no means confined to relativistic objects but are com- 
mon in the accretion discs that occur around young stars. In 
fact, not long after accretion disks were applied to quasars 
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and mini- quasars ( Lynden-Bcll, 1969, 197l|), and the galac- 



tic c entr e (Ekers & Lynden-Bell, 1971; Lynden-Bell fc Rees 
1972), hey w ere applied to star formation (Lynden-Bell & 



Pringle 1974). The jets that emerge from star-forming disks 
do not travel at relativistic speeds but at speeds of one or 
two hundred km/s. There is a correlation between the max- 
imum circular rotational velocity within the object and the 
spe ed of the jet and we give an explanation of this in sections 
2.3 



and 3.3 



First numerical calculations of the force-free 



structures we are advocating, were made by Li et al. (2001 ) 
but to date they have only calculated the first two turns 
whi ch give indications of the in itial growth of magnetic tow- 
ers. Krasnopolskij ct al. (1999|) studied the launching o f jets 



loaded with plas ma, as did Romanova et al. (1996 ), Cton 
topoulc s (1995a ) has - like us - looked for force-free self- 
similar solutions for jets and has considered t he possibility 
of purely toroidal field s (Contopoulos, 1995b ), Winds and 
Jets are considered by Lovelace et al. (1991) a nd Poynting 
jets from loaded winds have been computed by Lovelace et 
al. (2O0|2l). 



2 MAGNETO-STATIC THEOREMS AND 
DEDUCTIONS 

2.1 Theorems 

We consider force-free magnetic configurations with the feet 
of the flux tubes anchored in the accretion disk at z — and 
the magnetic configuration confined by an external coronal 
gas pressure p(z) which may depend on height z. Where 
there is magnetic field, there is no gas pressure. Earlier ver- 
sions of these theorems were proved only for the constant 
pressure case in paper II. 

The work done against the pressure in making a mag- 
netic cavity whose cross-sectional area at height z is A(z) 
is 



p(z)A(z)dz 



(1) 



This may be thought of as a potential energy. The larger the 
cavity occupied by magnetic field the greater is the pV en- 
ergy stored, an asset that could be drawn upon in recession 
by contraction. The total potential energy of the configura- 
tion stored in both magnetic energy and W p is 



W 



-j- 1 I (Bl + B 2 y + B 2 Z ) dV + W P 



(2) 



= W X + Wy + W Z + Wp. 



Contributions to these energies from slices at different 
heights z are given by 



B x dxdy 



over the area A, etc and 
w p = p(z) A(z) . 

Evidently, W x = J w x dz, etc, and W p = J w p dz. The min- 
ima of the potential energy give stable equilibrium config- 
urations. At any equilibrium (stable or not) the energy is 
stationary, so the work done by any small displacement 
consistent with the constraints is zero. We consider a small 




Figure 2. A slice of a force— free field structure. The twists about 
the axis are not shown. The theorem is true without assumptions 
of axial symmetry. To prove theorem I we consider a virtual dis- 
placement in which the slice dz is expanded by a factor of fi. At 
equilibrium an infinitesimal expansion will do no work. 



vertical displacement caused by expanding the slice of the 
configuration between z and z + Sz (see Figure 2) so that it 
now occupies the region between z and z + /j.8z. The region 
below z is unchanged; the region above z + dz is lifted by 
(fi — l)Sz. The work done by the pressure at the edge of the 
slice is 

— l)Sz5z(-dA/dz)p. 

which is second order in Sz and does not concern us; but 
significant work is done on the pressure elsewhere because 
the area A(z') that was initially at z has been moved to 
z + (p, — 1) 5z, so the area at z' is now 

A -fa- l)5zdA/dz . 

Thus the change in W p due to the displacement is 



AW P = ((jl - 1)6 z 



p(z')(-dA/dz')dz . 

roc 

p{z)A(z)+ A(z')(dp/dz')dz' 



(fx-l)6z 



To get the second equality we integrated by parts and used 
the fact that A vanishes at sufficiently great heights. This 
expression does indeed yield pdV when the pressure is con- 
stant. 

Magnetic fluxes must be conserved during the displace- 
ment so 



B x —* A* 1 B X , By — *• ji 1 B y , B z 



B z 



within the expanded slice whose volume increases by a factor 
/x. Evidently 



8nAW x 



8ttAW z 



B x dxdyjj,Sz 



B x dxdySz 



(M- 1 - 1) 



B x dxdySz 



B z dxdySz 



Collecting all the AW terms dividing by (fj, — 
1) 6z and letting a — > 1 we have theorem I. 

/>oo 

- (w x +w y ) +w z +p(z)A(z) + / A(z')dp/dz'dz' = .(3) 
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If we integrate over all z (by parts for the last term), we 
obtain 

POO 

- (W x + Wy) + W z + W p + / zA(z)dp/dzdz = . (4) 

Jo 

The integrals in both Q and (||) are negative when p 
decreases with height but vanish when p is constant. If 
d\np/d\nz = — s, then the final integral in (Q) is — sW p . Our 
next theorem is obtained by taking our virtual displacement 
to be a uniform lateral expansion with no vertical shift. Here 
we can not apply the virtual work principle unchanged since 
one of the constraints is that the vertical magnetic flux is 
frozen at the anchor points all over the disk z — 0. Our lat- 
eral expansion would violate this. However, there is an ex- 
tension of the principle in which the constraints are replaced 
by the forces of constraint in the equilibrium configuration 
and the virtual work then contains a contribution from the 
work done against the forces of constraint. In our case these 
forces are those due to the disk that balance the equilibrium 
Maxwell stresses. 

-(47I-)- 1 [B R B Z -R + B^BA + \ (Bl - B\ - Bl) z] 

per unit area of disk. The work done on the disk in a radial 
displacement £ = fiR, is 



VW = (4tt) - V 



B r B z R 2 d<f)dR 



evaluated in z — 0. 

The change in W p is just (jj 2 — 1) J pAdz since 
A — > fi 2 A, similarly B x — > fi~ 1 B x , B y — ► fi^ 1 B y and 
B z -f ^ 2 B Z . Since dV -> /i 2 dV we find 

AW = (fi~ 2 -l)W z + (/x 2 - l) W p 

with NO contribution from W x or W y since they are un- 
changed. The condition of equilibrium is no longer AW = 
but rather that any decrease in W must be balanced by the 
work done against the forces of constraint i.e., 

- [dAW/dtA^ = W . 

Thus we deduce Theorem II 



W Z ~W V 



iW 



(5) 



Many scientists brought up on the pinch effect find it strange 
that the toroidal component of the magnetic field has NO 
effect on the change of magnetic energy in such a global 
lateral expansion. In fact, there is no general tendency for 
the whole magnetic structure to shrink towards an axis. In 
axial symmetry Wx + W y = Wr + W<f> and none of them 
changes in a uniform contraction transverse to the axis. 
We return to explain the relationship of this result to the 
pinch effect in section 2.1. 

If, in place of the lateral expansion at all heights, we 
freeze the configuration below some height z and laterally 
expand everything above that height, then it is the magnetic 
stresses across the plane at height z that become our new 
Wo so we call their contribution Wo(z). Likewise, the W z 
and W p involved are integrated over the region above height 
z only. We shall call these quantities W z (z) and W p (z). We 
may then generalise (^) to give the exact result ([]) 



To localise theorem II consider a uniform lateral expan- 
sion that varies slowly with z so that R — > [1 + fJ*(z)] R~ We 
shall take fi to be zero in z = and to climb to a max- 
imum at zi before declining again to zero. Provided jj,(z) 
varies slowly enough the shear fi(z) will not produce signif- 
icant radial field from the displacement of vertical field. If 
we neglect this effect, the scalings are as for theorem II and 
we get a version of it localised near z\ (which we then re- 
place by z). Then our theorem III reads for z/R S> 1 and 
dB/dz < \B\/R 

w z ~ Wp (7) 

Wo does not contribute as there is no displacement on z — 0. 
(^) may also be derived by differentiating the result (^) with 
respect to z. That derivation demonstrates that (Q) is valid 
provided d/dz (^Wo(z)), is small compared with w z , which 
is so provided Br <C B z at height z. Unlike theorems I and 
II, theorem III is approximate and only true well away from 
the disk. If we integrated it down to the disk it would conflict 
with theorem II except in the special case when Wo is zero. 



2.2 Winding Makes Tall Towers 

We shall use cylindrical polar coordinates but we shall not 
assume axial symmetry. Then theorem I may be written 

Wr + = W z + (I - s)Wp 

adding W z + W p and using (theorem II) 



W = Wr + Wtf, + W z + W p 



(4-s)W p + Wo 

(4- s)W z + i(s-2)W 



Remember that for the constant pressure case s = 0. 
Here we shall assume (— dlnp/dlnz) = s < 4. We shall show 
elsewhere that when (— dlnp/dlnz) > 4 the magnetic field 
balloons off to infinity as in the highly wound pressureless 
case of paper I We assume that even after several turns the 
field structure resists being wound still further and that the 
work done per turn is asymptotically a constant. Further- 
more, we shall assume that the boundary term Wo which 
only involves the field on z — tends to a limiting value. 
Then in each turn W will be raised by a finite AW and 
AW = (4 - s)AW p = (4 - s)AW z . Hence W p and W z 
must increase without limit as the winding continues. Now 
W p = J p(z)A(z)dz can increase by increasing A(z) at given 
z or by increasing the height to which the whole configura- 
tion reaches. However, increase of A(z) for given z does not 
increase W z since B z is of order F/A(z) where F is the 
poloidal flux. Thus J B 2 dxdy tx F 2 /A(z). Since W p and W z 
have to increase together we deduce that the height of the 
whole structure increases for each turn of the flux anchor 
points on the disk. This argument is reinforced by the very 
crude estimates of field structure that follow and by the more 
accurate but more specific field models calculated later. 
When the configurations are tall, certain simplifications oc- 
cur in the theorems. Defining {B 2 } — A^ 1 J J B 2 dxdy the- 
orem III becomes 



(B 2 ) = 8tt P (z) 



(8) 



W z (z) - Wp(z) = iWo(z) = — 



B R B z R 2 d(j)dR . (6) 



where the average is taken at height z. Secondly, there is 
only a finite poloidal flux and as the system gets taller and 
taller the radial flux through a cylinder gets more and more 
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spread out. Thus Br becomes very small compared with B^, 
and B z . If Wr is neglected in theorem I, we have 



W 4, = w z + w p + J A(z')dp/dz'dz' . (9) 
In the pressure constant case we see from (^|) and ([]) 



lo^ = 2w z = 2w p 
so 

(Bl) = 2(B 2 ) = 16tv P . 
When p varies we define 



(10) 



a 2 = w„ 



A(z') (-dp/dz) dz 



then we have 

(Bl) = (2 - a) (Bl) = (2 - a) 8np 



(11) 



When the vertical Maxwell stresses (8tt) _1 (fif - B 2 - B 2 R ) 
are integrated over a cross section we find a net tension of 
w z — w^ — wr which becomes a force driving longitudinal 
expansion when that quantity is negative as in (g), ^ or 

©■ 



2.3 Crude Estimates Give Essential 
Understanding 

Let the total poloidal magnetic field emerging from an ac- 
cretion disk on z — be F. We shall assume that it returns 
to the disk at some larger radius so that it is anchored at 
both ends. Suppose that a typical field line reaches a height 
Z and that the tall magnetic structure has a radius R. Now 
each turn of the poloidal flux generates an equal toroidal 
flux which must pass through the area RZ. Hence after iV 
turns of the feet relative to one another the typical B^ is of 
order NF/RZ and the volume is tvR 2 Z so 

8ttW-V = N 2 F\/Z . 

The vertical flux F goes once up and once down so if it goes 
up in the inner R/\/2 and down outside that \B Z \ ~ 2F/nR 2 
and have 

8irW z =4F 2 Z/(tyR 2 ) 

Finally the flux passes through a cylinder of radius R/y/2 
radially so 



Br = F/ (V2ttRZ) , 
and hence 

8-kWr = F 2 / (2-kZ) . 

For a constant pressure p we have 

W p = pixR 2 Z . 

Hence 



W = — [(N 2 tv 2 + i) /Z + (4iT 2 + 8tt 3 P F- 2 R 2 ) Z] .(12) 

Minimising W over all values of R 2 we find 

ttR 2 = F(2np)i , (13) 



so the area is determined by F and p whatever the minimis- 
ing Z may be and the two terms that make up the coefficient 
of Z are equal at equilibrium. Minimising W over Z we find 



1 



z = Nn (r/Vs) 

with R already given by ( |l3| ) independent of N, this clearly 
shows that the height Z increases linearly with N ^> 1. 
Indeed, if we write N = (2n)~ 1 Qt where fi is the relative 
angular velocity of the flux feet, then 



4^2 



QRt 



(14) 



so the height of a steadily wound magnetic structure will 
grow with a velocity oc ili?. Notice that R is typically larger 
than the R of the inner foot of the flux tube at whose ra- 
dius fi is determined. Thus even with this excessively crude 
model we can see why the velocities of growing magnetic 
towers are directly related to the velocities in the accretion 
disk. 

In this section we have assumed thus far that our grow- 
ing towers do not have hollow cores with no field in them; 
consider however the top of a tower on axis. There the mag- 
netic field must splay out radially before descending down 
the exterior of the tower. Normally such division of one field 
line only occurs at a neutral point but here the magnetic 
field must resist the external pressure p so B 2 /8ty cannot be 
zero. Thus on the axis there has to be a most interesting 
cusp point. The cusp points downward and has a vanishing 
opening angle. The axial field line comes to the cusp at finite 
field strength and there splays very gradually at first, but 
then more rapidly with the field strength B /8n balancing 
the 'external' pressure which has invaded the axis and its 
neighbourhood, down to the level of the cusp. For points 
above such a cusp the tower will be hollow. Of course, it 
may be that this cusp is a mere dimple in the top of th e 



tower, however the MHD calculations of Li et al. (2001), 
who were only able to calculate the towers for the first two 
turns, suggest that the towers may be hollow over a signif- 
icant fraction of their height. Let us recalculate the energy 
assuming that our tower has a hollow core radius of R c and 
a maximum radius R m . Following our previous calculations 

B^ = NF/ [(R m - Rc) Z] , 

SnWtf, = N 2 F 2 n [(R m + R c ) / (R m - R c )] /Z . 

\B Z \=2F/ [n(R 2 m -R 2 )] , 

8TTW z = 'iF 2 Z/[jY(R 2 m -R 2 c )] . 

B R = Ff [2^ [i (R 2 m + R'i)]? z} , 

8-kWr = [F 2 / (2tvZ)] (R 2 m - R 2 C ) I (R 2 m + R 2 C ) , 
8-kW p = 8ty 2 P (R 2 m -R 2 C )Z . 

Writing x = R c /Rm we deduce 

2 \ n 

z- 1 - 



W = Fix' 



+ 



(1-x) 

4(l-x 2 ) 



Rm + 



(15) 



+8Tv 3 pF' 2 R 2 n (l - x 2 



Minimising over R m keeping x and Z fixed we find 



6 D Lynden-Bell 



^ = (2 7 r 3 p)-* J F(l-x 2 )- 1 

so at these minima we have for all x and Z 



Wo = -FV 



(1-1) 



, (i-« 2 ) 



+8 (2tt 3 P )' F _1 Z 



Minimising over all a; in the range < x < 1 we find that the 
quantity in square brackets increases throughout the range 
so the least value occurs at x = R c j ' R m = 0. Thus, at the 
level of this crude approximation, the towers are not hollow 
but filled with magnetic field. 



2.4 Stability 

Potential energy minima are always stable but if the min- 
imisation is carried out under the restriction of axial sym- 
metry there may exist asymmetrical distortions that lead to 
lower energy configurations. A useful insight was taught me 
by Moreno Insertis. Magnetic towers in tension are stable 
to sideways bowing. Like Euler struts, tall towers in com- 
pression are unstable. Such ideas require modification in the 
presence of an external pressure. The tower of fluid giving 
the ambient pressure is not unstable although it is under 
compression. Thus a magnetic tower that merely gives am- 
bient pressure will not be unstable this way, rather it is 
pressures above ambient that cause such instabilities. The 
modified stability criterion is w$ — w z < w p . This criterion 
is only marginally satisfied by our constant pressure config- 
urations for which the equality holds, so those will be pretty 
floppy, but when p decreases with height the buoyancy term 
provides stability. 



2.5 The Pinch Effect a Pressure Amplifier 

By showing that the energy in the B^ and Br components 
of the field did not change in a uniform lateral expansion we 
demonstrated that there is NO tendency for a highly wound 
magnetic structure to contract overall laterally. Where then 
does the Pinch Effect come from? To elucidate this consider 
a specified toroidal magnetic flux F$ contained between in- 
ner and outer cylinders of radii Ri and Ro both of height 
Z. Minimising the total energy over all possible B$(R) that 
give the flux, the minimum occurs with oc R -1 and the 
energy is then 

W* = iF$Z- 1 /{]n(R /R i )] . 

In uniform lateral contractions or expansions Ro /Ri remains 
constant, but if we fix Ro then W$ can be decreased by 
making Ri smaller. 

It is this that gives the pinch effect. Since B 2 ^ oc R~ 2 
the toroidal magnetic field acts as a pressure amplifier de- 
livering on the cylinder at Ri a pressure Rq/R? times the 
pressure it exerts on the outer cylinder. The whole pinch 
effect fails if there is nothing to push on at the outside 
and indeed the field would expand outward to larger radii 
reducing the pressure on Ri to zero. The Pinch Effect is 
a pressure amplifier; amplifying nothing gives nothing. 
In the presence of a B z the amplifier has less gain. By 
our theorems (87r) _1 J B z 2nRdR = n {R 2 po ~ RiPi) and if 
(87T)- 1 J B 2 z 2nRdR = r/nRfoo then pi = (1 - r])p Rl/R 2 so 



the amplifier only amplifies the excess pressure unbalanced 
by Bl 



3 CONFINED FORCE-FREE MAGNETIC 
FIELDS 

3.1 Basic Equations 

Inside the tower the magnetic fields dominate over any other 
forces so the magnetic field takes up a configuration that 
delivers no body force inside the jet, i.e., a force-free con- 
figuration with j x B = 0. Thus j is parallel to B and we 
may write j = aB where a is a scalar function of position. 
Now both j and B have no divergence so we deduce that 
B ■ Vq = which implies that a is constant along each line 
of force. We consider axially symmetrical systems in cylin- 
drical polar coordinates (R, <fi, z) in terms of the poloidal 
flux function P(R, z) which gives the flux through a circle 
radius R at height z. By flux conservation 



B z = (2-KRY 1 dP/dR , 
and 

B R = -(2irR)~ 1 dP/dz . 
We write 
B$ = {2nR)- 1 /3 , 
and then deduce 

B = VP x V(0/2tt) + (3X7(</>/2tt) , 



(16) 
(17) 
(18) 
(19) 



(3 is an axially symmetrical scalar function of position. Eval- 
uating the curl of (|f^) we find 

47tj = V x B = 

- [Rd/dR (R^dP/dR) + d 2 P/dz 2 ] V (<f>/2ir) + 

+V/3 x V(I>/2it . (20) 

The force-free condition j = aB now gives, cf both 

V/? = 4vraVP , (21) 

and 

Rd/dRiR-^dP/dR) + d 2 P/dz 2 = -4tt5/? . (22) 

From (§f|) it follows that the normals to the surfaces of con- 
stant P are the normals to the surfaces of constant P so f3 
is a function of P and 

/3' (P) = 4ttS . (23) 

Inserting this expression for a into ( p2| ) we get the basic 
equation for force-free equilibrium (we write Q(P) for (5(3') 

RdfdR{Rr 1 dP/dR) + d 2 P/dz 2 = -(30 = -Q(P) . (24) 

From @ we see that B • VP = so that P is constant 
along a line of force. Indeed the equations for the lines of 
force follow from the condition ds||B which gives 

dR Rdcj> dz 
Br Btf, B z 

or using the above expressions for B 

dR _ Rd(j> _ dz . . 

~ dP/dz ~ P(P) ~ dP/dR ' ( ' 

We shall use this later to work out the twists of field lines. 
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3.2 Field Structure Above a Differentially 
Rotating Disk 

Consider the tube of flux that rises within an inner circle of 
radius Ri(P) on an accretion disk and returns to it at some 
outer radius Ro(P). There will be a differential twisting due 
to the fact that Qi > fio and both are radius dependent. 
We define fi(P) = fli — do so f2(P) is the rate of differ- 
ential twisting of the field lines labelled P. At any time t 
when the twist angle has accumulated to be $(P) = Q(P)t, 
the field lines labelled P will rise to some maximum height 
Z{P) before heading down to re- intersect the disk at Rq(P). 
We are concerned to know ho w this total twist <& is dis- 
tributed over the height z. In Lynden-Bell (2001) I found 



1// = A?/A + (1-Ai)7(l-A) 



(29) 



that great progress could be made by simply adopting the 
idea that each field line had a certain twist per unit height 
<&(P)/Z(P), however I did not use the full power of the vari- 
ational principle to derive the equations. Equations ( p5| ) for 
the lines of force show that 



d<t>/dz = f3(P)/ (RdP/dR) 



(26) 



whereas P and /3(P) are constant along each line, RdP/dR 
is not. Now as the magnetic tower grows, its cross-section 
at any height will have some dimensionless profile / with 



P = P m (z)f(\) 



(27) 



Here P m (z) is the maximum value P achieves on a cross 
section of height z. Now Z (P m (z)) — z, so P m (z) is the 
inverse function of Z(P m ). A is the fractional area of the 
cross section at height z, A = ttR 2 /A(z) = R 2 /R 2 n where 
Rm(z) is the radius of the magnetic cavity at height z and 
A(z) is its area. From these definitions the maximum value 
of / is one and since P is zero on axis^J and at the tower's 
surface / will be zero at A = and 1. If / is independent 
of height the tower is said to be self-similar, however, in 
general the form of the function /(A) may depend on height 
although such changes may be slow except at the tower tops 
and near their feet. From ( p7| ) 

RdP/dR = P m (z)2A/'(A) = 2P/(dlnA/dln/) . 

Thus from @ 

d(j>/dz = iP" 1 /3(P)(dlnA/dln/) . 

Within any small height interval dz the twist of the 
line for force labelled P will have two contributions, one 
^P~ 1 P(P)(d\n\i/dlnf)dz as the flux rises through z at 
Ai and a second — ^P _1 /3(dlnAo/dln/)d,z as line descends 
through z at Ao- Thus the total contribution to the twist of 
the line from this height interval will be 



d$ = IP-^P) 



"din (Ai/Ao)" 


*-(£ 


d<f> 


„) 


din/ 




i dz 





dz, (28) 



as both Ai and Ao correspond to the same P and are at the 
same height, z, the value of P/P m (z) = / is the same for 
each. 

Thus we may regard Ai and Ao as the roots for A of the 
equation /(A) = /. A simple example will illustrate this. 
Suppose / is given by 



1 If the tower were hollow P would be zero for a whole area near 
the axis. 



This / is clearly zero at A = and 1, we have chosen the 
coefficients to ensure / = 1 at the maximum at A = Ai . Now 



f 



(Aj - A)(Ai+A-2AAi) 



[A? + A(l-2Ax)r 
and at A = Ai 

r = -2/[Ai(l-Ai)] . 

Hence near the maximum of / we find it is well approximated 
by 

/=1-(A-A 1 ) 2 [A 1 (1-A 1 )]- 1 , 

which is actually exact for all A when Ai = A. The roots for 
A are Ao = Ai ± [Ai(l - Ai)]* (1 - /)' and 



din 



ln/ = (l-/r 



/[Ai(l-Ax)]2 
[Aj-Cl-AxXl-/)] 



(30) 



the curly bracket reduces to 1 when Ai = \. Near the maxi- 
mum 1— / <C 1 so the expression becomes (I — /) i (A^ 1 — l)i . 
Away from the maximum Ai ~ A?/ [l + (1 — Ai) 2 /] and 
Ao ~ 1 - (1 - Ai) 2 /(1 + Ai/) so we find 

dln (f) //dln -^ 



[l-tl-Ar) 2 /] ' + 
+ (l-A 1 ) 2 /[i + (i-Ai) 2 /]" 1 . 



For Ai = £ and / = j this gives i| i.e., almost l| but for 
/ = ~itis||~lj and for / = j it is l|. Except for the 
integrable singularity of (1 — /) ~i at the top of each field 
line we deduce that in this simple example dln(Ao/Ai)/dln/ 
does not vary greatly. The integrable singularity at the top 
of each field line, P = P m is not a peculiarity of the example 
chosen. Returning to the general case and looking near the 
maximum / = 1 at A = Ai the roots Ao and Ai must behave 
as Ai ± c\/l — f so quite generally dln(Ai/Ao)/dln/ will con- 
tain the (1 — /)~ 3 factor near / = 1 We shall show in paper 
IV that in one limit / = — eAlnA and in the opposite limit 
/ = 2 A or 2(1 — A) depending on whether A is less than or 
greater than i. In both these cases dln(A2/Ao)/dln/ is near 
1 far from the top of the field line and in the latter case it is 
(I - f/2)- 1 which varies from f to 2 over the whole range 
0> /> 1. 

Notwithstanding such variations we shall start with 
the rough approximation that twist is distributed uniformly 
with height so that its distribution function 



g(P,z)dz = dz/Z(P) 



(31) 



This integrates to 1 over the range z — to Z(P) and gives 
so simple a result that it may be used as a starting point 
of a refined treatment that gives greater weight near the 
tops of field lines see section 4. With the distribution ( |3l| ) 
we calculate the B^, flux between z and z + dz. Each turn 
around the axis of an element of poloidal flux dP generates 
equal toroidal flux so the twist $<?(P, z)dz, gives a toroidal 
flux (2-K)~ 1 &(P)g(P, z)dPdz in the element of height dz. 
Adding the contributions over all the P that reach height z, 
i.e., those with P > P m (z) we find that the total toroidal 
flux per unit dz is 
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B&Rm — I Bj,dR — 



(27vy 1 ^(P)g(P,z)dP . (32) 



We shall find it convenient to work in terms of average 
fields at a given height, so the left-hand-side we re-write 
as B$(z)R m . Our former (-B^) will be related to (Py>) by 
some structural constant related to the profile of B^ with A. 
We define the ratio J 2 by 



J 



< Bl > _ 4 J (3 2 {P)\- 1 dX 
(S ) 2 ~ [/ f3(P)\-id\] 2 



(33) 



Were the tower hollow for R < R c then multiply the right- 
hand-side of j33| ) by (R m — R c )/(R m + \B Z \). Similarly for 
the 2 components we define \B Z \ cf ©, by 



\B Z 



and 



I 2 = i 



A- 1 



2nR\B z \dR = 2P m /A 



*• < = 



< Bj > 
\Bf 



> 1 



(34) 



(35) 



If the profile / depends on height then I 2 and J 2 will in 
general depend on height too, but for the tall towers gen- 
erated by continual winding we may expect that the pro- 
file form, /, to settle down to a typical one except near 
the tower top and bottom. Thus I and J will not vary 
strongly and indeed for self-similar towers they are dimen- 
sionless constants. In paper IV we find for the pressure con- 
stance case I = 1.359, J = .799 and for the linear case 
I = 1.179, J = 1.098. 

Combining (|) with (§3) and (|j| and setting R c = 



p 

~A~ 



< Bl >2 



(2-kvY 



(ttR 2 



(36) 



in which all the symbols may be functions of 2. 

We use ( ptj ) to write R m in terms of P m and p. Using 
( p2| ) to evaluate B^ via (poj), (|^) gives us the average B 2 , 
at given z 



(Bl) 



J 2 



{2Tr)- 1 $(P)/Z(P)dP 



Rn 



(37) 



(B 2 ) is given in terms of P m and R m by (|3fj); while Br is 
so small that we could neglect B R altogether, it is not un- 
interesting to see how its inclusion could correct the result. 
We shall make a very rough estimate. The radial flux up to 
height 2 = Z through a cylinder of radius Ri = R m /\/2 
will be F — P m (z). This gives a mean radial field of Br = 
[F - P m (z)/] (V^nRmZ). We therefore write 



(Bl) = [F- 



P m (z)] 2 I (2n 2 R 2 m Z 2 ) 



Substituting the above estimates into 8nW we have SttW 



^Att^PIR^ 2 + 8n 2 p{z)R 2 m + 



8nW 



J^(2n)-^/ZdP 



+nj 2 



+ (2n)- 1 (F- P m fZ~ 



+ 



> dz 



(38) 



Here P m {z) and Rm(z) are functions of 2 to be varied, 
Z(Pm) is the inverse function of P m (z) while p(z) and $(P) 
are given fixed functions. /, J, and F we treat as fixed con- 
stants. Only the first two terms involve R^, varying it and 
demanding that If be a minimum for all such variations 



gives the two terms equal, but that equality merely repro- 
duces what is already contained in equation (^) in the form 
nRm — (27rp) ~ "3 P m l ■ Inserting this R^ into the first term 
and doubling the result as the second is equal to it, those 
terms of 8tyW now reduce to J 4I^/87vpP m dz. Later we find 
it useful to change the independent variable from 2 to P m . 
We therefore write dz — dZ / dP m dP m and on changing the 
limits appropriately we obtain 



47 V /87Tp(Z) (-dZ/dPm) PmdPrr 



We now introduce II(Z) = j 41 yj9>irp(z)dz and after an 
integration by parts in which Z and II vanish at one or 
other end point we have just U(Z)dP m in place of the 
first two terms of 8ttW. 

Using Z(Pm) as our variable in the remaining terms we 
find that we do not have to vary both the function P m (z) 
and its inverse function Z(P m ), but W still contains both 
Z{Pm) and Z(P). They are the same function (albeit of a 
different variable). The third term of 8irW in (B8I) simplifies 
greatly when integrated by parts 



($/Z) dP 



'dZ/dP m )dPm = 



$(Pn) / §/ZdPdP m . 

Jo 

Again, we used the fact that Z(F) is zero in the bound- 
ary terms. We reverse the order of integration using 

So So™ (• • ') dpdp ™ = So Jp (■ ■ dP m dP to obtain 



*/Z / $(P m )dP m dP 



Finally we exchange the dummy variables P m and P to write 
the third term of W in terms of Z(P m ) 



$(P TO ) \Z{P„ 



$(P)dPdP„, 



(39) 



Inserting all these simplifications into 8kW and using 
C(-Pm) = Z~ x as our variable we find 



8nW= / ndP m + 
'0 



dP m + 



+ (F-P m f (27T)" 1 (d(/dP m ) dPrr. 



(40) 



where II = Tl(Z) may equally be considered as a function of 



3.3 Solution of the Variational Equations 

Varying £(P m ) and demanding that W be & minimum is 
now easy since the second and third terms are linear in £ and 
<5il/<5£ = —4IZ 2 tJ8ttp(Z). Thus the variational equation for 
Z[Pm) reads 



4IZ 2 y/8^pJZ) = (2n)- 1 J 2 ${P rn ) J QdP+ir-^F-Prr. 



.(41) 
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Figure 3. Profiles of the Magnetic cavity in a constant pressure 
environment. The differential rotation of the feet of the field lines 
labelled P is Q(P) oc (P + P)~ 7 with P small and 7 = 1.5. The 
profile of the cavity is given at two times. 
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Figure 4. The 'central dipole' model gives a more rounded profile 
than the power law model. 4a is at constant pressure, 4b has 
pressure falling as (1 + z) — 2 and is drawn for three times. 



Thus, the function Z(P m ) is given by 
Z[p(Z)]i = 



(tt/2); 



(4tt) 



o(p m ) // n{P)dPt 2 + 

+tv- 1 (F-P) 



.(42) 



The final F — P term arises from our rough estimate of 
the Br term; unlike the t 2 term next door, it does not grow 
with time, so after a few turns it is negligible as we expected. 
Neglecting it we define 



H(P„ 



( Pm ) Pjj 



tt(P)dP 



and we obtain the very pleasing result 

z\p(z)]i = CiPkn{p m )t , 



(43) 



(44) 



where 

Ci = (jr/2)* [J/(47r)]J-i . 

At given P,„ ( fjif ) gives Z as a function of t. Knowing that 
function we know how behaves with t. That then shows 
us how Z/t behaves with t] (^) shows us that, at constant 
Pm, Rm also behaves as Since p(z) is a given function 

this serves to define Z(P m ) or for that matter P m (z). Of 
course we still have to derive the values of the dimensionless 
structure functions I and J that we gave earlier, but these 



to 



only affect the time-scale of the evolution. Using 
re-express P m in terms of R m ( fli] ) takes the pleasing form 



t (J/7)i? m ni 



(45) 



which is a more sophisticated version of our result ([14|) . No- 
tice that (^) holds at each value of P m thus the collima- 
tion Z(P m )/R(P m ) grows linearly with time even when the 
pressure depends on height. 

The shapes of the magnetic cavities as functions of time 
follow directly by substituting the value of P m in terms 

of Rm and p(Z) from @ viz P m = [2^p(Z)\ * R'i/I 
into Cl(P m ) in ([ii]). After that J42] ) becomes an equation 
connecting the radius R m of the magnetic cavity to the 
height Z at which the cavity has that radius. We draw 
some examples of the evolution of the cavity shapes as 
functions of time in Figures 3 and 4. When p is inde- 
pendent of z, Rm(Pm) is independent of time by (Q) so 
all the expansion takes place along z axis. However when 

1 _ ! 

p decreases with height Rm(Pm) oc P m [p{Z)\ 1 and the 
last factor increases as z increases so it does increase with 
time. In spite of this Z/R grows linearly with t as seen in 
(Eq). A simple model will illustrate this, we take p(z) = 
po(l + z I a)~ s . We need s < 4 so that the field does not splay 
out as in paper I. Equation ([44] ) for Z now reads, writing 



V = (2tt)i [J/(27r)]P^n(P„ 
as a simple example 



= Vt. Taking s = 2 
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z = vt 



vt 



+ ^/l + V 2 t 2 /a 2 



Thus, when the pressure falls asymptotically like z~ 2 , Z(t) 
starts out at constant velocity V but then accelerates and 
asymptotically it behaves as 2V 2 t 2 /a of course ram pressure 
and inertia which are missing from our model may reduce 
this acceleration and ram-pressure could even reverse it, nev- 
ertheless, without such our model predicts acceleration when 
the pressure falls with height. The radius corresponding to 
P m is now given by 



Pm OC Pfr 



Vt 



+ y/l + VH 2 /c 



For small Vt/a this is proportional to 1 + Vt/a but for 
large Vt/a this factor changes to 2Vt/a so asymptotically 
Rm is proportional to time when Z is proportional to t 2 . 
These results are dependent on our choice of s = 2. For gen- 
eral s, Z asymptotes to Vt [Vt/a] 3 ^ 4 ^^ while R m asymp- 
totes to [Vt/a}] 3/(4 - a) . Both these expressions demonstrate 
the very rapid expansion as n approaches 1 i.e., when p 
falls almost as fast as z~ 4 and indeed our former consid- 
erations were limited to pressure that fell less fast than 
r~ 4 in spherical coordinates. The great extent of radio jets 
might however be an indication that we should be consider- 
ing not merely enclosed fields but those configurations that 
"reach out to infinity" We shall consider these again in an- 
other paper, but we remark here that if a field carrying the 
poloidal flux F reaches out to a region beyond which the 
pressure behaves as Hr~ 4 then the field will escape to infin- 
ity with B 2 — 8tyHt~ 4 with an opening angle theta given by 
2tt [1 - cos(0/2)] = F/V^rH. For small F/VH such open- 
ing angles can be small 

9 = [2F 2 /(tt 3 H)]^ 

and in such configurations the field will be close to radial 
from the source. It is possible that we should be considering 
the large scale radio jets as examples of this phenomenon 
rather than taking the pressure to fall off less strongly than 
r -4 , but currently I feel the other case to be more general. 



3.4 Estimation of the Twist Function P(P). 

In the discussion following equation ( j3o| ) we showed that 
the concentration of twist toward the top of each field-line 
behaved as (1 — /) _2 . If we substitute this behaviour for 
dln(Ai/Ao)dln/ in equation ( p^ ) and integrate along each 
line we find 

nF 

$(P) oc ^P _1 /3(P) J [1 - (P/P m )]"3 (-dZ/dP m )dP m . 
With Z(Pm) now known we can deduce the form of /3(p) 



/3(P) oc 2P$(P) 




-dZ/dP m )dP„ 



P (Pm-Pp 



(46) 



We can not go further without specific models for f2(P) and 
p{z). 



3.5 Displaced Power Law Models 

fi(P) has to be zero at P = F and it must fall as P increases 
from to F. If Q is finite at P = then the simplest 'power 
law' model is the doubly displaced one 

^^^[(Po + P)- 1 ~(P + F)-'] . (47) 

We calculate in the regime Po <C F, F — P ^> Po in which 
circumstance the (Po + P)~ 7 term is negligible. 

n(p ra ) r wdp = t^&^" 7 

Jp m (7-1) L 

-(Po + Pr' 7 " 1 '] , (48) 

where From @ the above expression 

is proportional to Z[p(Z)]~s . Two cases give power laws; 
7 > 1 gives Zp~z oc vj~^ J ~ i't, while < 7 < 1 gives 
Zp? oc ■cu~ 1 / 2 . We take our former displaced power law 
for the pressure with s < 4 

p = po [(z/a) + l]~ s . 

Notice that, for z < a/s, the pressure is almost constant 
so that regime is equivalent to taking s = but, when z is 
large, p oc z~ s so Z [p(Z)]~ oc Z 1 ' 4-3 ^ 4 . Hence we have the 
power laws 

Z oc zu t , 



(49) 



where 



r = 



(7 - i)S for 7 > 1 and S = l/[(s/4)] 



^jS for 7 < 1 

Rewriting (^) in terms of vj setting to* = Po + P 



p oc 2Pn(P)t 



(zu - P )t(-dZ/dtzr) 



(uj — VJ„)^ 

Now write x = zo/zu* and recall O(P) oc zu^ 1 ■ 

P oc Pzul - 7 Mf 1_s , 

where 

x- {P /zu,)]i x~ (T+1) dx 



dzu 



(50) 



(x- 1)4 



and we have assumed [(F + Po)/(P + Po)] 3> 1 to allow the 
upper limit to tend to infinity. 

Notice that, for P ^> Po, M is independent of P and 
indeed M = 7r 2 r!/(r - 1)! which is 3vr/4 for V = § and 
tends to 1 as F -> 1; for P = 0, M = 1. 

For the constant pressure case s = 0, S = 1, and when 
7 > 1 



/3ocP/(P + P)3 , 

so when Po is small compared with P, P oc pi 
j(S-l)-iS 7 >1 



In general T — 7 = 



For s = 2, S = 2 so F - 7 = 



17(5-2) 0<7<1 
'7-I 7 > 



(51) 



(52) 







< 7 < 1 



In paper IV we shall solve for the fields in such power 
law cases. Particular interest centres around cases in which P 
is proportional to P' or P which give rise to linear partial 
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differential equations for P c.f. equation (^). The former 
occurs when pressure is constant i.e., S = 1 when 7 > 1 
wherever Po is small compared with P. This has solution / = 
-eAlnA. The latter (J3 oc P) occurs when S = 27/(27 - 1) 
and 7 > 1 i.e., when s — 2/7 and with S = s — 2 with 7 < 1. 
Although the other power laws give non-linear equations 
they too can be solved except near the top of the tower. 
Again we leave the discussion to paper IV. Rotation in a 
gravity-dominated disk cannot fall off more rapidly than 
fid oc R~z while if the B z varies as R~ v on the disk then 
P oc R 2 ~ v . Thus fi d (P) oc p-s/C 4 - 2 "). The requirement of 
a finite field on axis suggests that near there r\ = and 
fid oc P~4 i.e., 77 = I but if the field further out falls off as 
R~ i then 77 = § and 7 = 3 so 7 is likely to be > 1 because 
B z is likely to fall faster than R~i 



3.6 Models with a Central Dipole 

If a central object of radius R a rotates uniformly at angular 
velocity fi s and carries a dipole of strength D then initially 
the field at radius R on the disk is B z — D/R? for R > R s . 
If the accretion disk is dominated by the central object then 
in centrifugal force gravity balance 

fi d = (GMs/R 3 )^ 

Now before it is twisted the flux function P of the dipole 
is P = 27rPr" 1 sin 2 so on the disk P = 2tyDR~ 1 for 
R > R s - Now suppose fi s < (GM S /R 3 s )i then there will 
be a co-rotation point in the disk at some point R — R c = 



where P = P c = 2-kDRZ 



[fi t/(GM a 

Then for P > P c we shall have a differential winding 
of the flux feet with fi s — fid negative while for P > P c we 
shall have a differential winding with fi s — fid positive. 

So 



fi(P) = fi s - (GM S )5 



(2ttD) 



fi s 



1 



(0 



.(53) 



Assuming that all the flux is wound up there will no doubt 
be continual reconnection close to P = P c where the toroidal 
flux runs in opposite directions. It then becomes of interest 
to ask whether a greater poloidal flux crosses the equator in 
the range R 3 to R c or in the range Po to 00 . Since P falls 
like P" 1 outside the star half the flux is within R q / W = 2R S 
so if the star rotates at less than \ // 8(GM B /Rs)i then the 
majority of the poloidal flux is dragged forward by the disk 
while if the star rotates at closer to breakup speed the disk 
will be dragging the majority of the flux backward. Of course 
the central object's field may be so strong that it enforces 
co-rotation as the inner disk without being able to do so 
further out. There is yet another interesting radius in the 
problem, if we think of the oppositely wound fields then at 
what radius is the rate of generation of toroidal flux within 
balanced by the rate of generation of net toroidal flux of the 
other sign on the outside? A little thought shows this can 
only occur when P c = (|)3P S = 1.84P a . On integration 
(||) gives, using Q with P c = F 



where y = P m /F 



J 2 / _1 Pyfi 2 



fi = fi 3 



{y 1 l-v-f } 



4 AN IMPROVEMENT IN ACCURACY 

We showed earlier that the twist of the field line labelled P 
is given by 

d^/dz = i/3(P)P _1 dln(A i /A ) /din/ , 

and that this last factor had an integrable infinity that be- 
haved as (1 - /) 2 = [1 - (P/P m )]"i. Here we show how 
to take some account of such a distribution of twist with 
height. In the total twist we shall have the factor 

z 

[l-{P/P m )T*dz = Z+ 



" (P/Prr. 



|(-dZ/dlnP m )dlnP„ 



The quantity within { } is small when (P/P m ) is small but 
becomes large close to P m = P. Near there —dZ{P m )/d\nP m 
can be approximated as constant at the value — dZ/dlnP. 
Since most of the contribution to the integral comes from 
near there we may take —dZ/d\nP m constant at that value 
and move it outside the integration which then becomes 
writing x 



£{P) = 



Pm/P and 1 

F/P 

'(1- 



V 



I) 



1 dx = 



V 1 - p /- F 



(l + y)' 1 2dy = 21n 1 + P/F 



The small region where P is close to F inhabits very lit- 
tle of the volume as these lines of force only just emerge 
from the disk before returning to it and, for P small, this 
integral is close to ln4 = 1.38. We may take account of the 
increased weight near the top of each field line by taking the 
distribution of twist with height to be the sum of a uniform 
distribution and a 8 function just below the top. Thus 

g(P,z)dz = [1 + S(z - Z-)(-dZ/dlnP)C(P)] N^ 1 Z~ 1 dz , (54) 

where Z- is Z(P) minus a very small quantity so that inte- 
grating up to Z(P) integrates over the 5 function. Ni is the 
normalising factor that ensures that g integrates to 1. 

Ni = [1+ (-dhiZ/dlnP)C(P)] . 

Now Z decreases as P increases so —d\nZ/d\nP is positive 
and will show no great variation. Since jC(P) is around unity 
for much of the range and Ni is bounded below by 1 we 
shall simplify by replacing — dlnZ / dlnP by —d\nZo/d\nP in 
the expression for Ni , where Zq is the approximate solution 
given by (H) . We insert the above g into (^2|) and find in 
place of (j33|) 



{Bl) = [J/(27t)] 2 



§P m C{P m )N- l Z- l + 
J Pm */{NiZ)dP 



R n 



Writing £ for Z 1 our expression for 8irW$ becomes 



(An)- 1 J 2 



$ 2 p^£ 2 ivr 2 dlt+ 

+2$P m £Af 1 |§- J Pm $(N^ 1 dP+ 



dP„. 
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The fii jst and last terms are linear in £ so their varia- (1970| ). Such ideas may have application to some 7-ray 



tion is simple. The central term while homogeneous of de- 
gree 1 in £ is not linear because of the dln£/dP m term. 
Varying this central term inside the integral and writing 
* = d/dP{$PC/Ni) we find 



-25C, { +Z<f(P m ) r o Pm */(JViZ)dP- 

-^"Vp 1n[Z(P)/Z(P m )]V(P)dP 



dP„. 



The coefficient of 8C, is homogeneous of degree zero in Z 
but causes significant difficulties when one tries to solve the 
equations for the minimising Z(P rn ). We surmount these 
difficulties by evaluating these terms with our zero order so- 
lution Z = Zo(Pm) as given by (|44[). If greater accuracy is 
required we can then solve and put back the solution itera- 
tively. 

The variation of SnW^ as given above leads us to the 
following equation for S(8ttW) in place of ( ffl| ) 

r-F 



4I^8ttZ 4 P (Z) 
= {27t)- 1 J 



N-^(p m )j; m ^ d p- 

-$iV 1 - 1 £P ra *(P m )- 
-$ 2 N^ 2 £P m ~ 



+ 



xln 



Zq(P) 



Z Q {P m ) 



dP 



(55) 



All terms on the right are known so this gives us a re- 
fined solution for Z{P m ). Every variable outside an integral 
is a function of P m and, except where stated explicitly, vari- 
ables inside an integral are functions of P. While the above 
expression may be more accurate it provides far less under- 
standing than our former solution (ff4|). As our aim should 
be insight rather than accuracy we have used ( fli) ) in pref- 
erence to dHB]). 

A method of improving on the solution ( [H^ ) so that the 
variable twist with height is properly allowed for is given in 
the Appendix. 



5 CONCLUSIONS 

A sequence of static models can be more illuminating than 
detailed dynamical simulations. The statics should be un- 
derstood before the dynamics is attempted. The continued 
winding of the magnetic field of an accretion disk will build 
up tall towers which make magnetic cavities, towering bub- 
bles of magnetic field in the surrounding medium. Such a 
conclusion does not depend on axial symmetry. 

Non-axysymmetri c behaviour is observ ed in the pretty 
plasma experiments of 



Hsu fc Bellan (2002) 



When that symmetry is imposed we can calculate the 
tall tower shapes of the magnetic cavities for any prescribed 
winding angles $(P) = £l(P)t and for any prescribed pres- 
sure distribution p{z). Examples of these towers are shown 
in figures 3 and 4. The heights of these towers grow at a 
velocity closely related to the maximum circular velocity in 
the accretion disk. Our primary results are encapsulated in 
equations (^) and (ff|). 

As stellar-mass black-holes form, the winding of the 
magnetic field of the collapsing core ma y cause jets to emerge 
from supernovae as first predicted by Le Blanc fc Wilson 



bursts, to the micro-quasars (Mirabel & Rodriguez, 199!:) 
and to SS433 flMargon, 198^ ). The possibility that the elon- 
gated hour-glass planetary nebulae, some of the most deli- 
cate objects in the sky, may be magnetic towers arising from 
the accretion disks of their central binaries is particularly 
appealing. 
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APPENDIX 

APPENDIX A: PERTURBATION THEORY 
FOR THE VARIABLE TWIST PROBLEM 

We would like to give the field lines the correct distribution 
of twist with height approximately 



g(P, z) = [l-(P/P m )]-*dz/L(P) , 
where 



(Al) 



L{P) 



[1- (P/P m (z))]-* dz 



[1 - (P/P m )] - 3 (-dZ/P m )dP„ 



putting this into (|35| ) gives in place of (^) 

' l x [1 -(&)]"* L-dP j 

which in turn gives an 8irW$ 

8ttW^ = (Any 1 J 2 x 



p-2 



Here the function to be varied Z(P m ) is deeply embedded 
in the integrals that define L{P m ) as well as occurring ex- 
plicitly. There is no simple way of solving the variational 
equations that result, that is why we chose to consider the 
simpler forms of g given in equations ( |3~ij ) and (^). How- 
ever the variational principle ( j4C| ) based on (^) is so simple 
that a perturbation theory based on it can be developed. We 
write W for the energy based on (A.2) and W for the energy 
based on fl39) so that 



W[t} = W{{] + (W, - W*) 



(A3) 



Now 8nW minimises at Co and only the first n term is non- 
linear so we may write 

8ttW[(] = 8ttW[( ] + 

+ [ ^ 2 n/sc 2 (C - Co) 2 + ^ 3 n/5< 3 (C - { ) a dP m . (A4) 

Jo 

Here Co = l/Zo which is the solution given by (^4TJ . We 
believe that Ci = 1/Zi where Z\ is given by mm will be 
close to the minimum of W. Our problem is to find the 
correction of Ci whic h will give 8W/ 8( = more accurately. 
Evidently from (A3) 



5W/5C, = 8W/8C, + <5(W* - W^/SC, 



Hence using (A.4) and evaluating 5(W$ — W$) etc., at the 
approximate solution we find 

C - Ci = - 

[8k5 {w* - w*) /sci + s 2 n/sci (Ci - Co) + ^n/scKd - Co) 2 

where we have written 8 2 H/8£q for 8 2 Tl/8C, 2 evaluated at Co 
etc. Now 811/ 'SC = -4IZ 2 y / 8-K P {Z) so 

8 2 n/8( 2 = AIZ 2 d/dZ \z 2 ^8np(Z) 



and 

8 3 n/8( 3 = -4IZ 2 d/dZ 



Z 2 d/dZ{Z 2 ^/8Tvp(Z)) 



so all those quantities are quite simple functions of Z and 
we have to evaluate them at Z = Zq the solution (^). Thus 
all terms on the right are known except 8nS(W,f, — W^)/8C,i 
which we have to evaluate at £i . We refrain from giving the 
gory details of the evaluation but 



= _ 27r f EZ 

where E(P, 



dP m + 8C AFmil 



J Pm ($/2ttL) [1 - (P/P„ 
d 



'5 dP and 



8Y/8C, 



HP) 



2ttL2[1-(P/P to )] 2 




(A2) 



In these terms we must put Z = Z\ the solution given in 
([55|). Thus we can in principle obtain by perturbation theory 
a solution to the variable twist problem. 

This paper has been typeset from a TfX/ LTjhX file prepared 
by the author. 



